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We study magnetic domain patterns in ferromagnetic thin films by numerical simulations for 
a simple Ising-like model. Magnetic domain patterns after quench demonstrate various types of 
patterns depending on the field sweep rate and parameters of the model. How the domain patterns 
are formed is shown with the use of the number of domains, the domain area, and domain-area 
(DX)' distributions, as well as snapshots of domain patterns. Considering the proper time scale of the 

& ' system, we propose a criterion for the structure of domain patterns. 
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I. INTRODUCTION 

in 

Various kinds of domain pattern are seen in a large number of physical and chemical systems and have been 
investigated experimentally, numerically, and theoretically 1,2 Magnetic domain patterns in uniaxial ferromagnetic 
garnet thin films also show many types of structure. For example, a hexagonal or square lattice, bubbles, stripes, 
and so on are observed under static or oscillating magnetic field perpendicular to the film plane. Under zero field, 

S after quench, a labyrinth structure is often observed. In some cases, the domain pattern under zero field has a sea- 
island structure, which consists of many small up-spin (down-spin) "islands" (domains) surrounded by the "sea" of 
down (up) spins. Labyrinth and sea- island structures are metastable patterns and should change to a parallel-stripe 
structure under the thermal fluctuation or other kinds of fluctuation. In fact, it was observed in experiments that a 
labyrinth structure changes to a parallel-stripe structure under field cycles;^ 
i ^ i • Magnetic domain patterns show different structures depending on the external field. There are several simulations 
which demonstrate how domain patterns change under very slowly changing field.— Some of them also showed 
hysteresis curves^ 7 . Though a hysteresis curve has no information about the structure of a domain pattern, it gives 
information about magnetic ordering. Magnetic ordering depends not only on the external field but also on the field 
changing rate. In fact, there are also some works about the dependence of hysteresis on the frequency and amplitude 
of an oscillating field£&^ 

\ In this paper, we focus on the domain patterns under zero field and study them by numerical simulations. Our 
model is a simple Ising-like model which reproduces domain patterns observed in experiments for ferromagnetic garnet 
£ — ' thin films. In experiments, domain patterns under zero field are observed after a strong applied field is switched off. 
Therefore, we need to start a simulation from the state under magnetic field above the saturation field in order to 
reproduce the domain pattern under zero field. Recently, the field sweep-rate dependence of domain patterns was 
studied experimentally and numerically in Ref. [ill However, the parameters except for the field sweep rate were 
y~\ fixed to specific values in the simulations. The parameters are different among samples in experiments. Here, we will 
also study a parameter dependence of domain patterns. The parameter dependence as well as the field sweep-rate 
dependence may give useful information on rich properties of the pattern formation in a ferromagnetic thin film. 

The rest of this paper is organized as follows. In Sec. HH we describe our model and numerical procedure. In 
Sees. IIIII and IIV1 the characteristics of domain patterns are displayed for fast- and slow-quench cases, respectively. We 
will show not only snapshots of domain patterns but also the number of domains, domain areas, and the domain-area 
distributions. Such quantities enable us to picture how the domain patterns appear. In Sec. El a criterion about 
^ ■ domain patterns and the field sweep rate is discussed. Conclusions are given in Sec. |VH 

_ _ i 

II. MODEL AND NUMERICAL PROCEDURE 



Our model is a simple two-dimensional Ising-like model whose Hamiltonian consists of four energy terms: uniaxial- 
anisotropy energy H anl , local ferromagnetic interactions (exchange interactions) Hj, long range dipolar interactions 
Hdi, and interactions with the external field g ex ^ i n i 12 i 13 We consider a scalar field <j>{r), where r = (x,y). The 
positive and negative values of <f>(r) correspond to up and down spins, respectively. The anisotropy energy prefers 
the values <fi(r) = ±1: 

fl „, = „/ dr (_«r)! + «rn. (1) 
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The exchange and dipolar interactions are described by 



Hj = p , ir i^pl (2) 



and 

tf d i = 7 J drdr'cf)(r)cj)(r')G(r,r'), (3) 

respectively. Here, G(r, r') ~ |r— r'|~ 3 at long distances. The exchange and dipolar interactions may be interpreted as 
short-range attractive and long-range repulsive interactions. Namely, <f>(r) tends to have the same value as neighbors 
because of the exchange interactions, and it also tends to have the opposite sign to the values in a region at a long 
distance because of the dipolar interactions. These two interactions play an important role in creating a domain 
pattern with a characteristic length. The term from the interactions with the external field is given by 

H ex = -h(t) J drct>(r), (4) 

where h(t) is the time-dependent external field. Here, let us introduce a disorder effect. In other words, spatial 
randomness is introduced in the model. We consider the disorder effect only in the anisotropy term: a is replaced by 
a\(r), where 

A(r) = l+/x(r)/4. (5) 

Here, /i(r) is an uncorrelated random number with a Gaussian distribution whose average and variance are and /tig, 
respectively. However, fi(r) should have a cutoff so that A(r) is always positive. In our simulations, we set /^o = 0.3. 
From Eqs. (U)-©, the dynamical equation of our model is described by 

90(r) 8{H &ni + Hj + H Al + H c ^) 

-L 



dt u 5<j){r) 

= L o |aA(r)[0(r) - <f>{r) 3 } +pV 2 cj)(r) - 7 J dr' '<f>(r')G(r,r') + h(t) 



(6) 



Hereafter, we fix Lq = 1. 

It is useful to calculate the time evolutions of Eq. ([6]) in Fourier space when we perform numerical simulations. The 
equation is rewritten as 

^ = a[(<f> - 3 )A] fc - (/3fc 2 + iG k )cjy k + h(t)S hfi , (7) 

where [•]& denotes the convolution sum and G k is the Fourier transform of G(r,0). If G(r,0) = l/|r| 3 , then 

G fc = a - aik, (8) 

where k = \k\ and 

a = 2tt / rdrG(r), oi = 2tt. (9) 

Jd 

Here, d is the cutoff length, which can be interpreted as the lower limit of the dipolar interactions. In fact, G(r,0) 
cannot be assumed to be l/|r| 3 at short distances, and Eq. ([8]) is the asymptotic form for k — > 0. However, the details 
of G(r, 0) is not important at short distances because the exchange interactions (the k 2 term) are dominant there. In 
the simulations below, we use Eqs. ((SJ and © as G k and set d = tt/2, which results in a,Q = 4. 

Equation Q is useful also for discussing some characteristic properties of domain patterns. Let us consider only 
linear terms in the equation 

^^=%0fc, (10) 
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where r/p. is the linear-growth rate for zero field. It means that <pk decays exponentially for negative r/k and grows 
exponentially for positive although the nonlinear term prevents </>fc's growing too much. The linear-growth rate 
has a quadratic maximum: 

Vk = -(Pk 2 - 7aifc + 7a ) + a 

/ \ 2 2 2 

= A 207 +^-7ao + «. (11) 

This suggests that the characteristic length of the domain patterns should be 2ir/ko, where ko = aij/2(3. Therefore, 
if the values of (3 and 7 are fixed, the characteristic length should be also fixed. In our simulations, we give those 
values as (3 = 2.0 and 7 = 2/3/ai = 2.0/tt, so that ko = 1.0. On the other hand, one of the central subjects in this 
paper is the a dependence of domain patterns. When a becomes large, the k region with positive r/k broadens and 
many modes of rjk grow. Therefore, the surface of domains can be rough for large a. 
In this paper, we consider a descending field as the external field 

h(t) = h ini - vt, (12) 

where hi n i is the initial external field, which should be equal to or larger than the saturating field. The value of the 
saturating field can be estimated from a simulation under ascending field as follows. After a domain pattern is formed 
spontaneously under zero field, an increasing field is applied. At a certain value of h, the domains with negative </>(r) 
disappear. Then, we consider this value as the saturating field. Although the saturating field depends on a, we set 
^ini = 1-5 for all values of a in our simulations, so that h ml is larger than the saturation fields for all a. 

The details of the numerical procedure is as follows. As the initial condition, <p(r) is given by random numbers 
in the interval 1.0 < 4>{r) < 1.1 at t = 0. The external field decreases from h- ln j to 0, following Eq. (fT2")) . Once the 
external field becomes zero at to = ^ini/w, the field remains zero for to < t < 2tQ. Namely, we stop the calculation 
at 2to- The time evolution is calculated by using a semi- implicit method. In other words, while we use the exact 
solutions for the linear and field terms, the second order Runge-Kutta method is applied to the calculation for the 
nonlinear term. The semi- implicit method enables us to use a rather large time interval: St — 0.1 in the simulations 
below. The simulations are performed on a 512 x 512 lattice with periodic boundary conditions. 
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FIG. 1: Domain patterns under zero field after a fast quench (v = 10~ 2 ), which display the dependence on the parameter a: 
[(a) and (b)] a = 1.5, (c) a — 2.0, (d) a — 2.5, (e) a = 3.0, and [(f) and (g)] a — 3.5. They are the patterns at t — 2to except 
for (a) and (f); these two are the ones at t = to- The size of the snapshots is 96 x 96, although the system size is 512 x 512. 
The value of (j>(r) is positive and negative in the white and black areas, respectively. 

III. FAST QUENCH 

In this section, we focus on the fast-quench case where the sweep rate of the field is v = 1CP 2 . The scale of the field 
sweep rate will be elucidated in Sec. [V] Considering the results in Ref. [ll|, we expect that some kinds of sea-island 
structure should appear in this case. It will be shown how the characteristics of domain patterns depend on the 
anisotropy parameter a. The characteristics are demonstrated by the snapshots of domain patterns, the total number 
of domains, average domain area, and domain-area distributions. 

Figure [T] shows the domain patterns under zero field. The white and black areas correspond to positive and negative 
(j>(v), i.e., up and down spins, respectively. For a = 1.5 and a = 3.5, the domain patterns considerably change after 
the external field becomes zero at t = to. By contrast, little additional change in domain patterns can be seen 
for a = 2.0, 2.5, and 3.0, so that the snapshots for them are displayed only for t = 2trj. When a = 1.5, black 
domains connect with each other due to the exchange interactions [see Figs. Ufa) andQJb)]. In the other cases (i.e., 
for a = 2.0, 2.5, 3.0, and 3.5), the domains do not tend to connect because of the dipolar interactions. However, 
the exchange interactions prefer connection of domains, and it is preferable in the viewpoint of energy. When a is 
small, the exchange interactions are relatively large and domains can connect with each other, though the coalescence 
of domains causes temporary energy loss of the dipolar interactions. For a — 3.5, after large black domains appear 
[Fig. [ljf)], their shape changes to be like an island with the same characteristic width as in other cases [compare 
Fig-Hfe) with Figs. EUcHUe)]. 

From Fig. [TJ one should notice the following: the larger a, the more inhomogeneous width of black domains. 
Namely, some domains are partly narrow or thick when a is large. The property can be explained by the linear- 
growth rate, Eq. (fTTj) . When a becomes large, the maximum value of 7]k also becomes large. Then, the k region where 
r/k > broadens. Therefore, ^fe's with different values of k grow and make the domain patterns with inhomogeneous 
thickness. The large black domains in Fig. [TJf ) are also explained in the same way. After some time, Fig.[ljf) changes 
into Fig.QJg) because a structure composed of domains whose characteristic length is 2ir/ko is more stable than the 
one composed of large round domains. 

The field dependence of the total number iVtot and average area A av of negative-^/ - ) (black) domains is shown in 
Figs. Ufa) and^b). Under a decreasing field, AT to t starts increasing at a certain field and stops increasing at another 
low field. In the plateau region of N tot , each domain grows and A m keeps increasing. When a = 1.5, AT tot starts to 
decrease after the plateau region. Moreover, for a = 1.5, N to t and ^4 av keep decreasing and increasing, respectively, 
after the external field becomes zero. This behavior is caused by connection of domains as mentioned above: cf. 
Figs.[lja) andQJb). If the calculation is continued for t > 2to, the values of TVtot and ^4 av will decrease and increase 
more, respectively. On the other hand, when a = 3.5, A av increases also after the external field reaches zero, although 
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FIG. 2: (Online color) The field dependence of (a) the total number of domains N to t and (b) the average domain area A av \ (c) 
the Aav-iVtot graph. The broken line in (c) expresses the curve of A av iVt t = 512 2 /2, where the total areas of black and white 
domains are equal. The calculation was performed under a rapidly decreasing field (v = 10 -2 ) until t — to and then under 
zero field until t = 2to- The number of domains means the number of negative-^(r) (black) domains. The domain area is the 
number of grid points that have negative values of c/>(r). 



Ntot does not change under zero field. This means that each domain just grows without increase of the number of 
domains. For t > 2to, A av does not increase when a = 3.5. In fact, A av of a — 3.5 has the maximum value at a 
certain time between to and 2to- 

Figure [2^ a) also enables us to know that the larger the a, the lower the value of the external field where the first 
black domain appears. This fact implies that spins cannot easily flop when a is large. In other words, for large a, the 
minima of -f/ani are low and the transition between negative- and positive-^(r) states does not easily occur. Another 
point we notice about Fig. [Ha) is that the larger the a, the smaller the N to t- This can be explained by the linear 
growth rate, Eq. (jTTJ) . As mentioned above, the k region where <pk > broadens when a becomes large. Therefore, 
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FIG. 3: (Online color) Domain-area distributions at t — to and t = 2to- The left, center, and right histograms at each time 
are for a — 1.5, a — 2.5, and a = 3.5, respectively. The domain area is the number of grid points that have negative values of 
<Kr). 



the size of each domain can become large for large a. If the size of each domain is large, the number of domains 
should be small. 

Figure [DJc) shows the A av -iVtot graph. The broken line is the curve of A av N to t = 512 2 /2, where the total numbers 
of the black and white domains are equal. The curve for each a starts at (A av , N to t) — (0,0) and approaches the 
broken line. The curve for a = 2.0 and a = 2.5 do not reach the broken line, which means that there is a remanent 
magnetization. 

Figure [3] displays the domain-area distributions at t = to and 2t for a = 1.5, 2.5, and 3.5. The distribution for 
a = 1.5 at t = io has a sharp peak at a certain small value of domain area. The peak falls and a long tail appears 
at t = 2to. This change reflects the connection of domains for a = 1.5 as mentioned above. For a — 2.5 or a = 3.5, 
the shape of the each distribution is similar to a Gaussian distribution and its change between t = to and t = 2io 
is small. Here, we note that the larger the a, the broader the distribution appears. It can also be explained by the 
linear-growth rate. 



IV. SLOW QUENCH 

In this section, we show our results for the slow-quench case where the field sweep rate is v — 10 -4 . Considering the 
results in Ref.llll we expect some kinds of labyrinth structure in this case. The dependence on a of the characteristics 
of domain patterns is demonstrated by the snapshots of domain patterns, the total number of domains, and average 
domain area. Here, domain-area distributions are not displayed since each domain can grow too long under the slowly 
descending field. 

Figure 3] shows the domain patterns under zero field. In contrast to the fast-quench case, the domain patterns for 
a — 1.5 do not change so much between t — to and t — 2to. It is because the connection of domains has almost 
finished by t — to- In fact, the connection of many small domains happens also in this case for a — 1.5 (see Fig. [5]). 
On the other hand, when a is large, the width of black domains is inhomogeneous. This property is the same as the 
fast-quench case. We should note another property: There are many branches of black domains for large a. It is 
related to the inhomogeneity of domain width. As black domains grow, they bifurcate at the place where the domain 
width becomes large. 

The field dependence of the total number A'tot and average area A av of black domains is shown on a semilogarithmic 
scale in Figs.[5][a) and[5{b). Also, in this slow-quench case, 7V to t decreases after the plateau region for a — 1.5 because 
of the connection of domains. Though iV to t decreases and A av increases after the external field becomes zero for 
a = 1.5, they change very little for t > 2to- Comparing Fig. O with Fig. [21 we find two things. First, the difference 
between both data for a = 1.5 is rather small. Second, the larger the a, the larger the difference. The dependence on 
the field sweep rate was discussed in Ref. |Tl| only in the specific case of a = 2.5. Although the discussion is basically 
applicable for other values of a, little difference appears for small a because spins can flop easily regardless of the 
field sweep rate. 

The inset of Fig. [5Jb) shows the field dependence of A av on normal scale. For a larger than 2.0, A av grows linearly 
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in the low field region. The linear-growth region of A av is within the plateau region of N to t- This means that the 
magnetization linearly decreases in that region as the external field decreases. 

Figure Eljc) shows the A av -iVtot graph on a log- log scale. The end of each curve is almost on the broken line which 
expresses the zero magnetization. In other words, the remanent magnetization is almost zero in this case. 



V. CRITERION FOR THE DOMAIN STRUCTURE AND STANDARD PARAMETER FOR THE 

FIELD SWEEP RATE 

Now, let us consider the dependence on the field sweep rate. Some discussion about it was given in Ref . [ill. Namely, 
sea-island and labyrinth structures appear under zero field for rapidly and slowly decreasing fields, respectively. Such 
behavior was explained by the concept of crystallization: For a fast-quench case, high "supersaturation" lowers the 
nucleation energy, while the nucleation energy is high for a slow-quench case^i When the nucleation energy is small, 
many domains appear at once and cannot grow so long. By contrast, when the nucleation energy is large, a small 
number of domains appear and can grow very long. However, these discussions were just qualitative ones, and no 
standard parameter responsible to the field sweep rate has been provided. 

In order to predict the domain pattern in experiments, a criterion about the dependence on the field sweep rate 
should be elucidated. Since the proper time scale of each sample is different, the balance between the proper time 
scale and the field sweep rate must be considered to give the criterion. The proper time scale is related to Lq of 
Eq. ©, while we set Lo = 1 in the simulations. First, we should scale out Lo and a parameter (3 in Eq. ©. If we 
write t = LoPt, then Eq. ([6]) becomes 



d(j>{r) 
di ' 

where a = a/ f3, 7 = 7 / /3, and 



a\(r)[4>(r) - 0(r) 3 ] + V 2 <^(r) - 7 / dr'<f)(r')G(r, r') + h(t), 



(13) 



7 /T\ h[ n [ v 



L0P 2 



(14) 



Here, let us call v = v/(Lo(3 2 ) a scaled sweep rate of the external field. Next, we have to decide the crossover sweep 
rate v c as well as the scaled one v c . Here, crossover means the crossover from a sea- island structure to a labyrinth 
structure. Then, the value of v/v c becomes a rough standard in distinguishing sea-island and labyrinth structures. 
If Lo and P of a sample could be measured in experiments, v c of the sample would be obtained from v = v/(Lq/3 2 ). 
Then, we can predict the pattern for an actual v from the value of v/v c . 

We need to consider how to distinguish sea-island and labyrinth structures. Let us introduce the average domain 
length as a quantity to distinguish them, which can be obtained by dividing the average domain area by the average 
domain width. For each a and v, the average domain width was calculated from the autocorrelation function at 
t = 2io- In Fig. [HI the average domain length is plotted against the field sweep rate for a = 2.0, 2.5, 3.0, and 3.5. 
Recalling the snapshots in Figs. [1] and QJ we can expect labyrinth and sea- island structures above and below the gray 




FIG. 4: Domain patterns under zero field after a slow quench (v = 10 4 ), which display the dependence on the parameter a: 
[(a) and (b)] a = 1.5, (c) a = 2.0, (d) a = 2.5, (e) a = 3.0, and (f) a = 3.5. They are the patterns at t — 2to except for (a) at 
t = t . 
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FIG. 5: (Online color) The field dependence of (a) the total number of domains iVtot and (b) the average domain area A av \ (c) 
the ^lav-A^tot graph. The calculation was performed under a slowly decreasing field (v = 10~ 4 ) until t — to and then under zero 
field until t = 2to. The data are displayed on a semilogarithmic scale in (a) and (b) and on a log-log scale in (c), while the inset 
of (b) shows the same data as its main panel on a normal scale. The broken line in (c) expresses the curve of AwAU = 512 2 /2. 

patched area, respectively. Here, we define v c = 10 -3 in this case, although it is not a clear crossover point. Then, 
the crossover sweep rate in another case is estimated by v c — v c Lq(3 2 — 2.5 x 10~ 4 x Lo(3 2 . 

Now, let us introduce a parameter k = log 10 (w/w c ). Then, we can expect that a sea-island structure appears for 
about k > 1 and labyrinth one for about k < — 1. In fact, the experimental results of Ref. [Til support the expectation: 
Sea-island structure appeared for 2 x 10 5 Oe/s and labyrinth structure for 10 Oe/s. Although we do not know the 
value of v c in that case, it is certain that the difference between the values of k for sea-island and labyrinth structures 
is more than 2 [= 1 — ( — 1)]. Therefore, we suggest that the parameter k should be a useful standard parameter to 
distinguish sea-island and labyrinth structures. 

Here, we explain a linear decay on the logarithmic scale in the field sweep-rate dependence of the average domain 
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FIG. 6: Dependence of the average domain length at t = 2to on the field sweep rate v. The average domain length is calculated 
from the average domain area and the average domain width (see text). The domain pattern looks like labyrinth and sea-island 
structures above and below the gray patched area, respectively. 




time (t) 

FIG. 7: Schematic picture about the characteristic times in the nucleation process. See text for the details. 

length (£ av ) in Fig. [6] Let us recall that the data in Fig. [6] are confined to the cases where the total number of 
domains does not change after the external field vanishes [see Figs.[2][a) and UK a)], ^av is inversely proportional to the 
number of domains (iV-tot), which is equal to that of the created nuclei (iV nuc i) since no connection of domains occurs 
in the cases investigated in Fig. [6] As N nuc \ is proportional to the degree of supersaturation (AF), the evaluation 
of £ av is eventually reduced to that of AF. Now, we consider two characteristic times needed to explain Fig. [6] the 
induction time (ri) and the preplateau time (T2). Figure [7] illustrates the important stages of nucleation, including 
those times. The induction time starts at ibegin, when the supersaturation begins, and ends at t\, when the first 
nucleus appears, i.e., T\ = t\ — ibegin- The preplateau time starts at t\ and ends at £2, when iV nuc i(= A^ot) begins 
to take the plateau value shown in Figs. Ufa), Ufa), and El i.e., r 2 = t% — 1\. Then, AF is divided in two parts, AFi 
and AF2: AFi and AF2 are the partial degrees of supersaturation during the induction time (ti) and the preplateau 
time (T2), respectively. 

We now proceed in estimating AFi by comparing the energies of states (A) and (B) without and with a nucleus, 
respectively: (A) all spins are up, i.e., 4>(r) = +1 everywhere; (B) there is a circular nucleus with radius Ro where 
spins are down among up spins. Namely, 4>(r) = —1 and +1 inside and outside of the nucleus in (B). In both (A) 
and (B), the anisotropy energy is the same: i? an i = — ctS/4, where S is the area of the system. In (A), Hj = 0, 
Hdi = (27r/d)7S l , and H cx = —hS. Therefore, the energy for (A) is 



H A = -aS/4 + (2ir/d)jS - hS. 



(15) 
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On the other hand, in (B), the dipolar interactions H^i can be estimated from the ratio of combinations of spin pairs. 
Since the ratio of combinations is given by the ratio of areas, 



H di = {2n/d) 1 S 



S-TrR 2 " 2 
S 



irRnV „S-ttR 2 ttR 2 



S 



- 2- 



S 



_ E>2 

= C27r/ ( /)-.S'- (>/,/)- .47ri2g(l-^ 



Then we have 



H B [t) = -aS/4 + Qw/dfrS - (27r/d) 7 • AirR 2 1 



s 



h(t)(S -2itR 2 ). 



(16) 



(17) 



Rigorously speaking, there is an additional surface energy in (B) due to the exchange interaction (Hj ~ 27r/3i?o), 
which is negligible in the present situation with Rq <C S. When the difference between a pair of energies Ha and Hb, 



AH x {t) =H A - H B (t) = (27r/d) 7 • AttR 2 q 1 - 



nR* 



-2irR 2 Q (h ini -vt), 



(18) 



is positive, it quantifies the degree of supersaturation. Noting the numerical evidence that the induction time depends 
on Lq and a but not on v, we consider n to be independent of v. Then, AF\ becomes 



where we used the equality AiJi (the 



AFi = AHx(tx) = 2ttR 2 t 1 v, 
0, i.e., there is no supersaturation at the beginning. 



On the other hand, the degree of supersaturation accumulating between t\ and t 2 is generally written by 



AF 2 



dt 



dt 



AH 2 (t) 



(19) 



(20) 



where AH 2 (t) stands for the energy difference between the present state and its ideal saturated one, which correspond 
to multinuclci versions of Ha and Hb in AH±, at arbitrary time during the preplateau time. We may consider 
AH 2 (t) as a monotonically decreasing function since the increase of nuclei reduces the degree of supersaturation and 
supersaturation should vanish at t = £2- Noting AHziti) = AFi and Ai?2(^) = 0, we have 



AFn = 



t2 dt (-±AH 2 (t) } =[A// 2 (/,) 



AH 2 (t 2 )} = AF 1 . 



(21) 



From Eqs. (JTSJ) and (|2T]|. AF = 2A_F\. Recalling the above mentioned relation (£ av oc iV nucl oc AF" 



we have 



log 



10 '-av 



C - log 10 v > 



(22) 



where C is a constant. Equation (f2"2")) shows the linear decay on the logarithmic scale in Fig. [S] While this scenario 
has still room for improvement, it provides an essential explanation of the crossover between sea-island and labyrinth 
structures. The improved scenario will appear somewhere in the future. 

Recently, a similar behavior to Fig. [5] was reported in quenched ferromagnetic Bose-Einstein condensations 
(BECs)J^ The number of spin vortices in the quenched BEC depends on the quench time. The following fact 
was demonstrated in Ref. Il4t For a slow quench, the spin state has nearly a single-domain structure and no spin 
vortices appear, and some spin vortices appear for a fast quench. In other words, the faster the quench, the smaller the 
domain structure. The behavior is really consistent with our results. It will be an interesting problem to investigate 
the common mechanism of the dependence on the quench time. 



VI. CONCLUSIONS 



We have discussed the magnetic domain patterns for several values of a in the cases of both fast quench and slow 
one. Except for a = 1.5, sea-island and labyrinth structures appear for fast- and slow-quench cases, respectively. 
When a = 1.5, domains connect with each other and form a labyrinth (stripelike) structure. On the other hand, 
when a is large, domains tend to have many branches and their width is inhomogeneous. We have also shown how 
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the characteristics of the domain patterns, i.e., the number of domains and the average domain area, change under 
decreasing field. Moreover, we have introduced the average domain length as one of the quantities which characterize 
domain patterns. The dependence of the average domain length shows that the change from a labyrinth structure to 
sea- island one is continuous against the change in the field sweep rate v. Despite this fact, with the use of the crossover 
sweep rate v c , we propose the following criterion: A sea-island structure tends to appear when k — log 10 (u/t> c ) is 
larger than about 1, and a labyrinth one does when k is smaller than about — 1. 
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